home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / cswap.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.1 KB  |  49 lines
  1.       SUBROUTINE CSWAP( N, CX, INCX, CY, INCY )
  2. *
  3. *     interchanges two vectors.
  4. *     jack dongarra, linpack, 3/11/78.
  5. *
  6. *     .. Scalar Arguments ..
  7.       INTEGER           INCX, INCY, N
  8. *     ..
  9. *     .. Array Arguments ..
  10.       COMPLEX           CX( 1 ), CY( 1 )
  11. *     ..
  12. *     .. Local Scalars ..
  13.       INTEGER           I, IX, IY
  14.       COMPLEX           CTEMP
  15. *     ..
  16. *     .. Executable Statements ..
  17. *
  18.       IF( N.LE.0 )
  19.      $   RETURN
  20.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  21.      $   GO TO 20
  22. *
  23. *       code for unequal increments or equal increments not equal
  24. *         to 1
  25. *
  26.       IX = 1
  27.       IY = 1
  28.       IF( INCX.LT.0 )
  29.      $   IX = ( -N+1 )*INCX + 1
  30.       IF( INCY.LT.0 )
  31.      $   IY = ( -N+1 )*INCY + 1
  32.       DO 10 I = 1, N
  33.          CTEMP = CX( IX )
  34.          CX( IX ) = CY( IY )
  35.          CY( IY ) = CTEMP
  36.          IX = IX + INCX
  37.          IY = IY + INCY
  38.    10 CONTINUE
  39.       RETURN
  40. *
  41. *       code for both increments equal to 1
  42.    20 DO 30 I = 1, N
  43.          CTEMP = CX( I )
  44.          CX( I ) = CY( I )
  45.          CY( I ) = CTEMP
  46.    30 CONTINUE
  47.       RETURN
  48.       END
  49.